Large-Deviation Properties of Sequence Alignment of Correlated Sequences.

The significance of alignment scores of optimally aligned DNA sequences can be estimated through the score distribution of pairs of random sequences. It is necessary to obtain statistics for the relevant high-scoring tail of the distribution. For local alignments of iid drawn sequences it has already been shown that the often assumed Gumbel distribution does not hold in the distribution tail, but has to be corrected by a Gaussian factor. Real DNA sequences were observed to show long-range correlations within sequences, which are not correctly modeled by iid random sequences. In this publication the large deviation method that was used in previous studies is applied to local and global alignment of such sequences with long-range correlations. We study the distributions over the full range of the support and obtained probabilities as low as [Formula: see text]. We show that again a correction to the Gumbel distribution is necessary to study the dependence of the parameters on the correlation strength. For global alignments the Gamma distribution, which was found heuristically to be a good fit in earlier simple sampling studies, is found to be a poor fit.

miR-96 is required for normal development of the auditory hindbrain.

The peripheral deafness gene Mir96 is expressed in both the cochlea and central auditory circuits. To investigate whether it plays a role in the auditory system beyond the cochlea, we characterized homozygous Dmdo/Dmdo mice with a point mutation in miR-96. Anatomical analysis demonstrated a significant decrease in volume of auditory nuclei in Dmdo/Dmdo mice. This decrease resulted from decreased cell size. Non-auditory structures in the brainstem of Dmdo/Dmdo mice or auditory nuclei of the congenital deaf Cldn14-/- mice revealed no such differences. Electrophysiological analysis in the medial nucleus of the trapezoid body (MNTB) showed that principal neurons fired preferentially multiple action potentials upon depolarization, in contrast to the single firing pattern prevalent in controls and Cldn14-/- mice. Immunohistochemistry identified significantly reduced expression of two predicted targets of the mutated miR-96, Kv1.6 and BK channel proteins, possibly contributing to the electrophysiological phenotype. Microscopic analysis of the Dmdo/Dmdo calyx of Held revealed a largely absent compartmentalized morphology, as judged by SV2-labeling. Furthermore, MNTB neurons from Dmdo/Dmdo mice displayed larger synaptic short-term depression, slower AMPA-receptor decay kinetics and a larger NMDA-receptor component, reflecting a less matured stage. Again, these synaptic differences were not present between controls and Cldn14-/- mice. Thus, deafness genes differentially affect the auditory brainstem. Furthermore, our study identifies miR-96 as an essential gene regulatory network element of the auditory system which is required for functional maturation in the peripheral and central auditory system alike.

Score distributions of gapped multiple sequence alignments down to the low-probability tail.

Assessing the significance of alignment scores of optimally aligned DNA or amino acid sequences can be achieved via the knowledge of the score distribution of random sequences. But this requires obtaining the distribution in the biologically relevant high-scoring region, where the probabilities are exponentially small. For gapless local alignments of infinitely long sequences this distribution is known analytically to follow a Gumbel distribution. Distributions for gapped local alignments and global alignments of finite lengths can only be obtained numerically. To obtain result for the small-probability region, specific statistical mechanics-based rare-event algorithms can be applied. In previous studies, this was achieved for pairwise alignments. They showed that, contrary to results from previous simple sampling studies, strong deviations from the Gumbel distribution occur in case of finite sequence lengths. Here we extend the studies to multiple sequence alignments with gaps, which are much more relevant for practical applications in molecular biology. We study the distributions of scores over a large range of the support, reaching probabilities as small as 10^{-160}, for global and local (sum-of-pair scores) multiple alignments. We find that even after suitable rescaling, eliminating the sequence-length dependence, the distributions for multiple alignment differ from the pairwise alignment case. Furthermore, we also show that the previously discussed Gaussian correction to the Gumbel distribution needs to be refined, also for the case of pairwise alignments.